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Abstract: We present updated results for anisotropy searches in the direction of the Galactic Center 
(GC) at energies in the ranges: 0.1 < E < 1 EeV and 1 < E < 10 EeV. We use data from the Pierre 
Auger Observatory up to March, 2007. Present analyzes are therefore based on a substantially larger data 
set than our previous published results. A limit on the flux coming from a hypothetical point-like neutron 
source at the GC for 1 < E < 10 EeV was imposed, and searches for extended excesses were also 
performed. 



Introduction 

The widely accepted idea that the Galactic Cen- 
tre harbors a super massive blackhole which would 
be associated to the radio emissions from Sagit- 
tarius A*, together with the recent measurements 
by the H.E.S.S. collaboration of TeV 7-ray emis- 
sions close to the location of Sagittarius A* (T), 
makes the GC a very interesting place to look 
for localized excesses of cosmic rays at EeV en- 
ergies (1 EeV = 10 18 eV), especially for Auger, 
where it passes only 6° away from the observa- 
tory zenith. Also, there have been claims in the 
past from the AGASA and SUGAR experi- 
ments of large excesses in this region, none of them 
were confirmed by a later analysis of Auger data 
J4J- Recently, several scenarios predicting neutron 
fluxes from the GC have been discussed in the lit- 
erature [5 1, so that it is interesting to update pre- 
vious Auger bounds with the larger accumulated 
statistics of the surface array. The bulk of cos- 
mic rays consists of charged particles, which are 
strongly deflected in the galactic magnetic field. 
We therefore perform as well an extended source 
search around the GC for windows sizes of 10° and 
20°. We divide our data set into 2 energy bands: 
0.1 < E < 1 EeV and 1 < E < 10 EeV. Both 
ranges are studied in terms of Li-Ma overdensities 
significances. 



Data Set 

For the present analysis, we use data collected by 
the Pierre Auger Observatory from January 1st, 
2004 until March 31st, 2007, representing a con- 
siderable enhancement in the statistics compared to 
previous Auger published results. We consider sur- 
face detector (SD) events where 3 or more tanks in 
a compact configuration are triggered. The events 
are required to pass the level 5 trigger J6), for 
which the station with the highest signal has to 
be surrounded by an hexagon of working tanks. 
Such a requirement ensures a good reconstruction 
of the event. Finally, only showers with zenith an- 
gle 6 < 60° are used. The energy assignment for 
each shower is done by correlating a SD energy 
estimator (in this case, a variable related to the 
time integrated water Cerenkov signal at 1000 m 
from the core) with the calorimetric energy mea- 
surement performed by the fluorescence detector. 
This is the procedure used to build the Auger spec- 
trum from SD data and it is discussed in details 
in 0. The final data sample energy distribution 
is plotted in figure Q] For the events used in this 
analysis, the angular resolution |8] is about 2.0° 
for 0.1 < E < 1 EeV, where the sample is dom- 
inated by showers triggering only 3 stations and 
1.2° for 1 < E < 10 EeV, with the enhancement 
of higher multiplicity events (see mean values of 
the histograms in figure |2j. 
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Figure 1 : Energy distribution of the events after se- 
lection cuts showing the low energy SD threshold. 
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Figure 2: Angular resolution (1.5er) distribution 
for Q.K E < 1 EeV and 1< E < 10 EeV. 



Searches around the Galactic Centre 

To impose upper limits on the flux of a hypotheti- 
cal neutral particles source at the GC, we have done 
the search using a Gaussian filter matching the ex- 
periment angular resolution, since this leads to an 
optimal search as shown in (9). We have assumed 
the following differential form for the CR spectrum 



acceptance as possible, like detector deadtimes, ar- 
ray growth and weather effects. Therefore, to actu- 
ally impose the limit, the acceptances A s (E) and 
Acr(E) for the source particles and for the CR 
bulk, as a function of energy, must be well under- 
stood, especially at low energies, where the SD de- 
tector is not fully efficient, as can be seen from the 
strong suppression close to the SD threshold (see 
figure Q]). 



/ E \ 



-3.3 



* CR {E) = K Wi — \ EeV-WV^sr- 1 0.1 < E < I EeV 



which has a slope consistent with the HiRes exper- 
iment in the energy range 10 175 < E < 10 18 5 eV 
and is a smooth extrapolation of the Auger spec- 
trum measured above 3 EeV. The factor k is in- 
troduced to take into account the limited knowl- 
edge of the true CR flux and should be of order 
unity for the existing measurements at EeV ener- 
gies (k — 1.2 if the HiRes normalization is the cor- 
rect one, around 2 for AGASA and 1 for Auger). 
The 95% CL upper bound in the integrated (and 
acceptance weighted) source flux is then given by 

Ernax E max 

[dEA s $f(E) = -^—Aira 2 ( cIEAcr&cr{E), 

E min E min 

where n® 5 is the signal counts upper bound at 95% 
CL for a given observed/expected ratio and n exp 
is the corresponding background counts inside the 
Gaussian window of size a. The background es- 
timation is done using coverage maps, which take 
into account as much of the local effects over the 



Below 1 EeV, neutrons from the GC are strongly 
depleted by beta decays. An interesting prospect 
for data below 1 EeV would then be the possi- 
bility of imposing upper limits on the flux of ul- 
tra high energy photons. The H.E.S.S collabo- 
ration has reported a remarkably fiat spectrum of 
gamma-rays above 165 GeV from Sagittarius A*. 
A naive extrapolation of this spectrum would lead 
to a flux of photons of 0.58 km~ 2 yr _1 and 0.14 
km -2 yr _1 above 10 17 eV and 10 175 eV, respec- 
tively. In table [T] we show the observed/expected 
ratio for a point-like search using a Gaussian beam 
with a — 1.3° and the corresponding 95% CL 
bound on the source counts. No significant excess 
is seen. As discussed above, the conversion of n 9 s 5 
into a photon flux upper bound needs a dedicated 
study of the photon acceptancf]. 



1. For the case of photons, not only we know the ac- 
ceptance to the bulk (nucleons and heavier nuclei) is dif- 
ferent from the source, but also the energy assignments. 
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search 


window size 




/ 95 
l^obs I Tlexp T^s 


extended 


10° (TH) 
20° (TH) 


5663/5657 = 
22274/22440 = 


1.00 ± 0.02(stat) ± O.Ol(syst) 
0.99 ± O.Ol(stat) ± O.Ol(syst) 


point-like 


1.3° (G) 


192.1/191.2 = 


1.00±0.07(stat)±0.01(syst) 17.9 



Table 1: Summary of excess searches below 1 EeV around the GC (Sagittarius A*) in the form of extended 
sources using top-hat (TH) beams and point-like sources with Gaussian (G) beams. A 1% systematic 
contribution from possible weather and detector deadtime induced rate modulation is included. The Gauss 
beam matches the angular resolution (AR=1.5er = 2° (8)) in this energy range. 408607 showers selected. 



search 


window size 






nf 


$f (km^yr- 1 ) 


extended 


10° (TH) 
20° (TH) 


1463/1365 = 
5559/5407 = 


1.07 ± 0.04(stat) ± O.Ol(syst) 
1.03 ± 0.02(stat) ± O.Ol(syst) 






point-like 


0.8° (G) 


16.9/17.0 = 


0.99 ± 0.17(stat) ± O.Ol(syst) 


5.6 


0.018k 



Table 2: Same as tableQ]for l < E < 10 EeV, including also 95% CL flux upper limits on a neutron source. 
The Gauss beam used matches the angular resolution (AR=1.5er = 1.2° [8 1) in this energy range. 109101 
showers selected. 



We have done as well a search for extended ex- 
cesses around the GC on top-hat windows of 10° 
and 20° and the results are presented in table Q] 
showing no significant excesses. The isotropy of 
the sky around the GC for 0.1 < E < 1 EeV can 
also be seen from figures |3]C) and D), which are, 
respectively, a map of overdensity significances 
and the histogram of such significances, together 
with the expectations of isotropy. Given the statis- 
tics at each bin, we expect that the Poisson distri- 
bution can be fairly approximated by a Gaussian, 
which is plotted in red. 

1< E < 10 EeV 

Given the neutron lifetime at rest (885.7s), a con- 
siderable fraction of neutrons produced in the 1 < 
E < 10 EeV energy range at the GC would sur- 
vive and reach the Earth. Even though energies 
above 3 EeV are probably too high to be reached 
by particles accelerated by galactic objects, in this 
energy range, the bulk of CRs will be below that 
value. Unlikely the case of photons, if we assume 
here that the composition of the bulk of the CRs in 
the energy range studied is similar to the source, 
so that the source spectrum $ s has the same spec- 
tral index as the CR spectrum ^ci?,, we can obtain 
directly the 95% CL upper bound in the integrated 
flux (see numbers in table [2] where we also show 



the results of extended source counting) as 

„95 m r a 
$95 = ^ 47r(7 2 / dE <S, CR ( E)i 

T^exp J 

which can be finally converted into an upper limit 
on the source luminosity of Cgc < 1.25k x 10 34 
erg/s. Such a limit already excludes most of the 
models for neutron production at the GC [5 ]. 

In figure [3] A) one shows the Li-Ma iflOl signif- 
icance distribution of overdensities in 5° top-hat 
windows in the region around the GC for the en- 
ergy range 1 < E < 10 EeV. The distribution 
of the significances of figure [3] A) is shown in fig- 
ure [3] B), together with the expectations from an 
isotropic sky. The overdensities are consistent with 
an isotropic sky. 

Conclusions 

With an exposure two times larger than the one 
used in the previous Auger published result, we 
have found no significant CRs flux excess in the 
region around the Galactic Center in 2 different en- 
ergy ranges: 0.1 < E < 1 EeV and 1 < E < 10 
EeV. We have treated the GC both as an extended 
and a point-like source. The distribution of Li- 
Ma overdensity significances in this region shows 
consistency with an isotropic sky for both energy 
ranges. For 1 < E < 10 EeV, an upper limit on the 
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Figure 3: A) and C): Maps of the CR overdensity significances near the GC region on top-hat windows of 
5° radius for 1 < E < 10 EeV and 0.1 < E < 1 EeV, respectively. The GC location is indicated with a 
cross, lying along the galactic plane (solid line). B) and D): Histograms of overdensity significances for the 
same energy ranges. The dots represent the data, the solid histogram is the MC prediction of an isotropic 
sky and the expected Gaussian fluctuations are given by the red solid curve. 



flux of a hypothetical point-like neutron source in 
the direction of Sagittarius A* was imposed, which 
excludes most of the models predicting neutrons 
from the GC. 
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